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ABSTRACT 





U 


A new method for predicting the optimal velocity distribution on the wall of a two-dimensional 
diffuser is presented in the paper. The method by Principle of Dynamic Programming solves the 
optimal control problem with inequality constraints of state variables. The physical model of 
optimization is to protect the separation of the boundary layer while getting to be maximum pressure 
ratio in a diffuser of a specified length (or getting to the shortest length in a specified pressure ratio). 
The calculation results arc fairly in agreement with the experimental ones. It shows that optimal veloci- 


ty distribution on a diffuser wall should be as: the flow decelerates first quickly and then smoothly, 
while the flow is near separation but always protects from it. The optimal velocity distribution can 


directly be used to design the contour of the diffuser. 


INTRODUCTION 

A diffuser is an important part of compressors, fans and other air ducts. More and more atten- 
tions have been paid to its design. In the past dozens of years, the popularization and development of 
the optimization technique make it possible to design a diffuser with optimal velocity distribution 

The index of optimizing a diffuser is to obtain the highest pressure ratio under the condition of a 
minimum constructional length. Generally speaking, in order to get an optimal shape of a diffuser, it is 
necessary to know an optimal velocity distribution on its wall. With the distribution, the boundary lay- 
er can be avoided seperation and a maximum pressure ratio (or pressure recovery) can be obtained in a 
specified length. 

Nowadays, most designs of diffusers, which are two-dimensional or axial-symmetrical, arc still 
based on experience. Designers often use the criteria of the diffusing angle or the equivalent diffusing 
angle and one-dimensional calculational method to design it. Obviously, it is too simple to reach the 
index of the optimal design. 

Stratford(1959) proposed that the loss in a diffuser with the minimum length is the minimum 
while the boundary layer inside it is close to but just before occurrence of separation, then the velocity 
distribution is the best and the shape of the diffuser is optimal. Some researchers, such as H. licbcck, 
H. Fcrnboly, have used this principle to make some optimal designs. 

Many authors also investigated the flow field in a diffuser and study how to control the flow 
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scparaction. 

Some authors attempt to use the optimal control theory to solve the optimal design of a diffuser, 
because the governing equations of the flow in it are differential ones. 

Gu and Ji (1987) proposed an optimal design problem of a diffuser, using the optimal control the- 
ory and the boundary layer theory. The optimal velocity distribution on its wall was obtained by using 
Pontryangin's maximum principle. 

In order to meet the demands of engineering application, the optimal problem has to satisfy some 
constraints in both aerodynamics and strength which can be divided into two parts: one is called as 
constraint of state variable and another constraint of control variable. Those constraints arc often 
inequality and make it very difficult to solve the problem in mathematical treatment. 

It is well known that Pontryagin / s principle can only solve the optimal control problem with con- 
straints of control variables. To overcome the difficulty, many authors have done some research work 
and modifications such as continuous transfer techiquc (Jocobson, 1969; Gu, 1987) and expanded pen- 
alty function method (Gu and Miao, 1987). 

However, for a problem with more inequality constraints of state variables the treatment is not ef- 
ficient which have been stated by Gu and Miao(1987). That is to say, The more the constraints, the 
more the difficulties. On the contrary, the principle of dynamic programming is quite good at treating 
of state and control constraints. The more the constrains, the faster the calculation, because the numcr 
of considered states and decisions decreases in seeking optimal decision. 

In the present work, a physical model and a mathematical expression for dynamic programming 
are established and calculated. The result yielded by the method is quite in agreement with not only the 
experimental ones but also the result by Pontryagin's maximum principle. 

ESTABLISHMENT OF AN OPTIMAL D ESIGN PROBLEM OF A DIFFUSER 

It is well known that the flow losses in a diffuser mainly consist of separation loss and friction one. 
Obviously, the former is greater than the latter. The friction loss is always inevitable. However, it 
doesn't vary greatly because the friction coefficient is approximately a constant in the fully-developed 
turbulent flow. The total friction loss can be considered as increasing proportionally with the axial 
length of a diffuser. So the key to designing an efficient diffuser is to avoid the separation of boundary 
layer. Considering the two factors mentioned above, Stratford. (1959) proposed that the properties of a 
diffuser with the minimum length and without boundary layer separation is optimal. 

Generally, the turbulence degree at the inlet of a diffuser in engineering is so high that we can as- 
sume for convenience that the boundary layer has become a turbulent one at the inlet edge. And 
incompressible flow is only considered in present work. 

The typical expressions of a optimal problem for a diffuser are as follows: 

A) Pressure rise coefficient is maximum (i.e. the discharge velocity is minimum) under the condi- 
tion of a provided constructional length and without separation of the boundary layer. 

B) The length of a diffuser is minimum under the condition of a provided discharge velocity (i.e. 
the pressure rise coefficient is known) and without separation of the boudnary layer. 

It can be proved that the expressions A) and B) are correlative. 
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The separation of the boundary layer is a very complex problem. According to the change of ve- 
locity in main flow, the separation of boundary layer can be predicted by some experimental formulas 
to some extent. The following equations are adopted as the basis of solving the optimal problem, (Ref. 
1,3,4,5,6,7,8,9,10). 

The velocity shape factor of boundary layer is introduced as follows: 


T = 


6 2 du 
Udx 



m = 1/4 
m = 1 


where 




turbulent flow 
laminar flow 


( 1 ) 


U is the velocity of main flow. 

r>0 denotes acceleration flow, and T<0 deceleration one. So r can be used to judge whether the 
separation happens or not. 

For a deceleration flow, the relationship between i5 2 and U is as follows (Rcf.l) : 



-0.0,75-4.15^ 
U dx 



1/4 


( 2 ) 


For convenience, the length of a diffuser, L, is used as a characteristic length, the velocity at the in- 
let, C, as a characteristic velocity. Then Eq. (2) can be rewritten in a non-dimensional form. The 
non-dimensional length of a diffuser contour is S = X / L, the non-dimensional velocity V = U / C, 
the non-dimensional momentum thickness 9 = S 2 / L. So the Re d 2 * s 


Re, = U 


< 5 2 / v = V 


0 • Re, 


where Re 0 = C • L / v at the inlet of the diffuser. 
Then Eq. (2) becomes as: 


0 dV 

= 0 . 0175 - 4.15 • -—(Re. 

V ds ‘ 


V • 9)' 


(3) 


(4) 


Substituting the non-dimensional form of Eq. (1) into Eq. (3), we yield 

— = (0.014 - 3.52 • D(Re • V - Q)~ U * 
ds 0 

According to the result of Nikuradsc's experiments, Buri proposed that the boundary layer will sepa- 
rate when r is not greater than — 0.06. As stated by Gu and Ji(l 987) to ensure the flow in a diffuser to 

be far from separation, we utilize the limit of T as : 

- 0.04 T < 0 


MATHEMATICAL EXPRESSION OF THE OPTIMAL PROBLEM 

The mathematical expression of the index of optimization A) is as follows: 

Index function: 

J(* ) = V(* )-» min 


s.t. 
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dV 3/4 Q-5/4 -1/4 

— = r • V *0 'Re 

ds 0 

~i 

(6-1) 

^ =(0.014- 3.52 • r)(Re •F-0) -1/ " 
ds 0 

(6-2) 

o>r^ - 0.04 

(7) 

K(o)=l, 0(o) = 0 o 

(8) 

^(*), Q(*),.free 

(9) 


So the mathematical expression for dynamic programming solution is as follows: 
Index function: 


/(l) 




ds -> min 


s.t. 


dV/ ds = Tv' / '§~ UA Re~ WA 


1/4 


dQ/ds = (0.014 - 3.52D (Re o VQ) 
initial condition: V(o)=l 0(o) = 0 o 

control constraint: — 0.04 ^ T < 0 

state constraint: 0<V<1 0 ^0 

o 

Quantizing the equations listed above, we yield: 

r(k)v(k) 3 /A m~ 5 /A Re~ WA 


J = mini X AS * 


n*) U-o 


V(k + 1) = V(k) + AS * r(/c)K(/t) 3/4 0(fc) W< Re o W * ] 

0(lt + 1) = 0(k) + AS * { (0.014 - 3.52r(*)[/?e o K(fe)0(Ar)l ~ W * } 


- 0.04 T(Jt) < 0 

0 < V(k) < 1 

0 4£0(/c) 

0 

and the iterative relation becomes 


J(F,0,fc)= min AS* [T(k)V(k) W *Q(k)~ iM Re o W *] + J(V,B,k + 1) 

rt*) ^ a 

J(V,Q,N) = Q 


CALCULATION RESULT AND ANALYSIS 


(10) 


( 6 - 1 ) 

( 6 - 2 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 
(19) 


( 20 ) 


In calculation, Reo and N are taken as 10® and 10 respectively. The state variables and control 
variables are quantized respectively. The sets of admissible state variables are as follows: 

V= (0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1.0} 

0— { 0.00226, 0.00508, 0.00781, 0.01058, 0.01336, 0.01890, 0.02168, 0.02445, 0.02722, 0.03} 

Then the allowed quantized states are: 
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{(0.4, 0.00226), (0.4, 0.00508), (0.4, 0.0078), "VVVV.O, 0.03)} 

And the set of admissible control variable is as follows: 

T = { -0.04, -0.035, -0.03, -0.025, -0.02, -0.01 5, -0.01 , -0.005 } 

The two-dimensional dynamic Programming Computatinal procedure is used because there are 
two state variables, the calculation procedure is presented in the computer flow chart. The initial value 
of V and 0, that is , the values at stage o, are taken as 1 and 0.00226 respectively. The calculation re- 
sults at stage 5 are presented in Table 1. There is only a part of all results because the results arc too 
many to list them all. At each stage, a similar table can also be listed. 

In the table 1, a grid point stands for a allowed quantized state, the value put to the right-up of a 
grid point is the optimal value of index J at this state, and the one put to the right-down of the grid 
point is the corresponding optimal value of control variable T. 
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As shown in Fig. 1, curve ABC is an optimal velocity distribution of the diffuser. The results are is 
quite in agreement with the experimental data and the results calculated by Pontryagin's maximum 
principle by Gu and Ji,(l 987). 

The optimal velocity distribution ABC is also called as an optimal deceleration curve for a 
diffuser. In the range above the curve, including the curve, there is no separation while T is not less 
than -0.04; on the contrary, in the range under the curve, separation will happen. So the curve ABC is 
a critical line (r = -0.04) between separation and non-separation. 

Drawing a deceleration curve AB'C' in the nonseparation range, we can find that for a diffuser 
with a specified length, the discharge velocity V * is always less then V'. On the other hand, for a speci- 
fied discharge velocity V' (which menas a specified pressue rise coefficient) the corresponding optimal 
length S* is always less than S'. That means reducing velocity along the optimal deceleration curve 
ABC will yield the maximum pressure rise coefficient in a specified length or the minimum length of a 
diffuser in a specified pressure rise. In this case, the length is minimum and the loss is nearly minimum 
because there is no separation. 

From the physical explanation of the optimal deceleration curve in Fig.l It is also proved that the 
two expressions of index functions A) and B) of the optinal control problem are correlative. That is, 
the optimal velocity distribution obtained by one index function can satisfy another automatically. 


THE CALCULATION OF THE DIFFUSER CONTOUR 


The contour of the diffuser is calculated by means of the optimal velocity distribution on the sur- 
face, so that it is also called as optimal design problem or inverse problem. It is well known that solving 
the problem directly in X-Y plane will involve non-identifed of calculated region. So coordinate 
transformation is neccessary. It is the easiest way to transfer the X-Y plane to 4>-'F plane. 

The governing Equations, in fl>-T plane has been deduced strictly in the paper as: 

Taking an element in X-Y plane and considering incompressible, potential flow, the continuity 


equation and non-rotation equation arc as follows: 
dV aV 

— + L = o 

ax 3y 

dV aV 

— y - =0 

ax ay 

Velocity vector is 

V = \ V\ • (i • cos/? + j • sin/?) 

Where |V| is the amplitude of V. /? is the angle between Kand coordinate line X. 
The transform relation between X-Y plane and <h-T plane is: 

f y = (-x v f 9 + x 9 -fJ/Ji 

Where J is Jacobi matrix. 

Substituting Eq. (23) and Eq. (24) into Eq.(21), because of J=£0, then we yield: 
cosj? — Y,,sin0 ) + %(-Y 9 cos0 + J*sin/? ) 


( 21 ) 

( 22 ) 

(23) 


( 24 ) 
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acos/J . 


asm/? 


- r. 


acos/? 
► aY 


) = 0 


( 25 ) 


The four items on the left side of Eq. (25) a re as follow s respectively, 
9 -^r • (Y^ • cos0 - A^sin^) = jj£ • + Y ^ 


dO 

dV_ 


<- r. 

( -x. 


asin /? 


ad> 

cos/? + A^sin/l = 0 

asinff y acos^^j = Q 
afl> 


* 90 ) 
3C0S^> ^ __ y 


•L . JF 

yi A Q 


+ r. 


® aT ' aT v -- ® ■ - ® 

In deducing the defmations of normal and tangent unit vectors of equal d> and equal *P lines, Eq. 
(25) can be rewritten as, 

(26) 




+ Y^ +V 


*§_ 

aT 




+ r, -o 

c> 


I t can be easi ly p roved that: 

-Jx‘,+ r’.-v/J 

And substituting it into Eq.(26), finally, the continuity equation in (P-T plane can be given as: 

alnK Bj _ (27) 

ad) + aT 

In the similar way, the non— rotation equation in d>— T plane is: 

alnF _ a£ (28) 

a¥ ad) 

Two Laplace's Eqs. can be obtained from Eq. (27) and Eq.(28): 

a 2 lnF a 2 ln F = Q (29) 

a Jv F ad> 2 

. Ml „ 0 (30) 

ad) 2 aT 2 

The velocity distribution within the diffuser can be obtained by solving Eq. (29) with ADI 
method. Then from Eq.(28) the values of fi on the top line (*P = 0 or *E = 1) of the potential flow region 
can also be yielded, so the shape of potential flow region can be defined. The diffuser contour can be 
modified by adding thickness of boundary layer. The calculation result is shown in Fig.4, the line of 
Y/ L = 0 is the central line of diffuser. The shape is quite similar to the real size of B. S. Stratford's dif- 
fuser. The equivalent diffusing angle of the diffuser is 19 ° (integral angle), and is much greater than 

ordinary recommended angle. 


CONCLUSIONS 

The optimal deceleration curve (i.e. optimal velocity distribution) on the wall of a diffuser is first 
obtained by using the principle of dynamic programming. In solving optimal control problem of fluid 
mcchnics with inequality constraints of state and control variables, the dynamic programming method 
has many advantages over others. The physical model of optimization for a diffuser is to avoid the 
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separation of boundary layer while getting to the maximum pressure rise in a diffuser of a specified 

length (or getting to the shortest length in a specified pressure rise). The calculation results are fairly in 

agreement with the experimental ones and the results calcaulated by Pontrayagin's maximum principle. 

The optimal velocity distribution on a diffuser wall should be as: the flow decelerates first quickly 
and then smoothly, and the flow is near separation but always protects from it. The optimal velocity 
distribution can also be expanded to design an unsymmetric diffuser. 
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